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ABSTRACT 

I discuss the use of cross correlations between galaxies with distance information 
and projected weak lensing dark matter maps to obtain a fully three dimensional dark 
matter map and power spectrum. On large scales I < 100 one expects the galaxies 
to be biased, but not stochastic. I show that this allows a simultaneous solution of 
the full 3-D evolving galaxy bias and the dark matter power spectrum simultaneously. 
Within the photo-z information of the CFH lensing legacy survey, this allows a three- 
fold reduction of statistical error, while a cross correlation with CLAR or other deep 
spectroscopic surveys allows a tenfold improvement in dark matter power accuracy 
on linear scales. This makes lensing surveys more sensitive to the cosmic equation of 
state and the neutrino masses. 

Key words: Cosmology-theory-simulation-observation: gravitational lensing, dark 
matter, large scale structure 



1 INTRODUCTION 

The great success of modern physical cosmology has been 
driven by precision measurements of the cosmic microwave 
background (CMB), whose properties can be computed to 
high accuracy from first principles. Current CMB fluctua- 
tions have been measured to high statistical accuracy, which 
describes the universe at recombination. To constrain funda- 
mental cosmological parameters, such as the matter density 
of the universe, one requires measurements at low redshift 
as well. 

The low redshift measurements are difficult. Galaxy 
distributions can be measured in three dimensions and to 
high accuracy, but their formation is not understood in any 
detail. This leaves a systematic error in using their power 
spectrum to infer cosmological parameters. While the CMB 
power spectrum has been measured to statistical accuracies 
better than 1%, galaxy power spectra are not reliable trac- 
ers of matter at any accuracy close to that. It is known 
that galaxies of different morphologies are biased relative to 
each other. To complicate matters further, it is now known 
that the bias depends on scale, and that their d istributio n is 
stoch astic on non-linear scales llTeemark fc et al 2003.: .FanI 
l2003fl . This bias and stochasticity is expected to evolve with 
redshift. The interpretation of their distribution is difficult 
to quantify, and galaxy power spectra can only constrain 
funda mental parameters to limited precision. The WMAP 
team llSoergel et alJl2003h combined the CMB power spec- 
trum with the 2dF galaxy power spectrum. For the latter. 



E-mail: pen@cita.utoronto.ca 



the power is measured to high statistical accuracy, but with 
an unknown bias, so its interpretation is limited by system- 
atic errors and was not used in their analysis. The WMAP 
team instead used only the slope of the galaxy power spec- 
trum, which is hoped to be a less biased indicator of the 
underlying dark matter distribution. But it is not indepen- 
dently known how large an error that might have, so the 
error on the error bar is very difficult to quantify. Using no 
external data, WMAP data allows even closed cosmological 
models with no cosmological constant. 

In contrast, the distribution of dark matter can be pre- 
dicted to high accuracy through perturbation theory and 
simulations, since their equations of motion are simple and 
can be modelled to high precision with N-body simulations. 
Its measurement is intrinsically very challenging. Statisti- 
cal gravitational lensing of background galaxies in blank 
fields has been a very successful approach. But the sig- 
nals are intrinsically weak, limiting the accuracy of the con- 
straints. Nevertheless, the cleanness of the data has provided 
some of the strongest low redshift cosmological constraints 
JVa,n Waerbeke et aUEonllContalHi et aflEua. 

One intrinsic limitation of gravitational lensing studies 
is the broad redshift distribution of lenses. The observed 
shear signal integrates over all the matter distribution along 
the line of sight. Using multiple source planes gives some 
handle on the distance, but even in the best of cases onl y 
very crude redshift resolution can be extracted|Huj^o3). 
Sample variance of lensing surveys, and fundamentally the 
cosmic variance of an all sky survey, are limited to the num- 
ber of two dimensional modes on the sky. If one could beat 
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this cosmic variance, one could significantly improve on the 
results. 

Galaxies, on the other hand, allow measurements of ac- 
curate redshifts. One can measure their redshifts from opti- 
cal or radio spectroscopy. Their sample variance is smaller 
by the volume to surface area ratio of a survey. But as men- 
tioned previously, the distribution of galaxies is known to be 
biased, and this bias is likely to depend on scale and cosmic 
time. Their evolution is unlikely to provide a direct me asure 
of the cosmic equation of state. IJain fc Tavlorl i2003ri have 
proposed a purely geometric procedure to measure the geom- 
etry of the universe. Our proposal differs by actually solving 
for the matter distribution, which allows much higher po- 
tential accuracy. 

In this paper I show how one can combine the redshift 
resolved distribution of galaxies to significantly improve the 
information of lensing surveys. It is possible to simultane- 
ously measure the evolution of bias as a function of scale 
and time, and thus the 3-D distribution of dark matter. The 
only intrinsic limit is stochasticity between galaxies and dark 
matter. 



2 STRATEGY 

Perhaps this sounds like a free lunch. Just counting degrees 
of freedom, the evolution of bias as a function of scale and 
time is a two dimensional function. The observable dark 
matter power spectrum is a one dimensional function of an- 
gular scale. The cross correlation, however, does have two 
dimensions as well. The map data set is also two dimen- 
sional on the sky, so the number of degrees of freedom adds 
up that one could in principle measure the dark matter evo- 
lution using galaxy cross correlations. 

To illustrate the procedure, we will consider a simpli- 
fied case where the three dimensional quantities are pro- 
jected along a Cartesian axis, say z. The analogy of this 
axis is distance along the line of sight, which also corre- 
sponds to the direction along which we expect cosmological 
evolution to be important. One generically expects the bias 
to depend on cosmic time. If bias did not depend on the 
length scale tangential to the line of sight, the procedure 
is straightforward. We assume that we know the three di- 
mensional distribution of galaxies up to an unknown bias 
function b{z), 5g{x,y,z) = b{z)5{x,y, z). Any spatial depen- 
dence of t he bias will be captured by the cross correlation 
l|Pen 1998) coefficient, which we will discuss later. The two 
dimensional distribution of dark matter is projected by a 
known lensing kernel so the lensing surface density 

K is K{x,y) = J S{x,y, z)'WL{z)dz. For this example, we ne- 
glect the change of angular scales with distance. We can 
form a galaxy projection weighted by an unknown weight 
function w{z). For practical purposes, this weight function 
will be smooth, so it suffices to sample it at a number of 
points riz . The galaxy surface density in terms of this yet to 
be determined function is 



T,g{x,y) = y 5g{x,y,Zi)b{zi)w{zi)^Zi 



(1) 



We define a residual function e = J(Sg — K)^dxdy. One 
can then solve for the weight function by difi^erentiating for 



weights at redshift Zi. This results in a coupled system of 
equations 



dxdySg(x,y,Zi)T,g{x,y) = / dxdyhi{x,y)Sg(x,y, Zi). (2) 



The unknowns WLizi) appear implicitly linearly in equa- 
tion @ through the surface density Eg and equation Q. 
One can solve this system of equations to obtain w{zi) = 
WL{zi)/b{zi) if the intrinsic distributions of galaxies and 
dark matter are the same, there is no noise, and the num- 
ber of pixels ripix in the map is larger than the number of 
unknowns in the redshift dimension. Typically, the bias is a 
slowly evolving function of time. One would bin in coarse in- 
tervals, for example a dozen photometric redshift bins. The 
image has many effective pixels. With the exact solution, 
the weighted galaxy map will be identical to the dark mat- 
ter map. 

In the presence of noise, we can write Equation @ as 
a linear system 



Aw = K zb o" 



(3) 



where A is an ripix by Uz matrix. We define a noise matrix 
N. If the primary source of noise is shot noise in the lensing 
map, this matrix is diagonal with entries being the variance 
at each pixel a^{x,y). One solves the minimum variance 
solution to ^ as 

w = (A*N"^A)"^A*N"^K. (4) 

In general one has ripix ^ n^, so the inverse operation is 
well defined. 

One can in a similar fashion compute scale dependent 
bias evolution. We define the scale dependent bias as the 
ratio of the power spectrum of galaxies to that of the dark 
matter. 



b^{k,z) 



-Pgal(fc,z) 
-Pdm(fc, z) ' 



(5) 



The power spectrum is a two point statistic, which depends 
on the square of the bias. In terms of the density distribu- 
tion, this relates the Fourier transform of the galaxy density 
field linearly through b(k, z) to that of the dark matter: 



5g{k,z) = b{\k\,z)5{k,z), 



(6) 



or equivalently states that the galaxy density field is a 
convolution over the dark matter distribution Sg{x,z) = 
J 5{x, z)b{\x\, z)d?'x. A bit of care needs to be taken in the in- 
terpretation of this Fourier transform. We have represented 
the three spatial coordinates {x, y, z) by a vector x. The ra- 
dial direction z is now the third component of the x. The 
three dimensional expression in equation Q assumes that 
b{z, r) varies sufficiently slowly as a function of z that one 
can apply a transform at each interval Zi . This is reasonably 
accurate assumption except for the largest scale, where one 
would want to solve a general linear parametrized system 
analogous to instead of Fourier transforming. 

The density can be thought of as a one-point distribu- 
tion, which depends linearly on the bias. For this redshift z 
and scale k dependent bias b{z, k) we can Fourier transform 
from wave number space k to separation r: 



biz, 



(2^)2 



, , , , sin(fcr-) , 2 ,, 
biz.k) — r — -k dk. 
kr 



(7) 
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Our redshift binned bias function Q now generalizes to 
a two dimensional function of redshift and separation. We 
then convolve each redshift slice by the parametrized b{z,r). 
We then solve for the spatial and temporal structure of the 
bias. In practice, the bias should evolve slowly in space and 
time. Once we know the evolution of bias as a function of 
length scale and cosmic time, it is straightforward to obtain 
the dark matter power spectrum as a function of scale and 
time. In the linear regime, this measur es the growth fac- 
tor. I n the non-linear regime, the PD96 jPeacock fc Doddj 
Il99(f) approach gives a mapping between linear and non- 
linear power spectra. This needs to be calibrated against 
simulations, but as we see below, the inferences are proba- 
bly systematically limited by other non-linear effects. The 
growth factor in turn maps directly onto the equation of 
state of sp ace-time a nd matter, i.e. dark energy and neu- 
trino mass teJ2002h . 



3 LIMITATIONS 

The improvement is intrinsically limited on several fr onts. 
It re quires that bias be linear and non-stochastic 
Il998fl . On non-linear scales , the power spec t rum of galaxies 
is known to be stochastic JPen et al.ll2003l: iHoekstra et alJ 
|20o3)- In the cross correlation tomography, one can mea- 
sure the degree of stochasticity which provides a fundamen- 
tal limit to the accuracy. Other sources of limitations are 
the noise in the lensing and galaxy power measurement. 

One ca n measure the impact o f stochasticity. We use the 
model from "Hock stra et al.1 liooj) . The power spectrum of 
galaxies is related to the dark matter by the bias, Pg = h^P, 
and the cross correlation coefficient r(fe, z) can depend on 
time and scale. It relates {SgS) = r^JPgP. At the two point 
level, the bias and cross correlation completely describe all 
two point statistics, and parametrizes all effects including 
non-linearity. 

In the presence of stochasticity, the solution using our 
procedure becomes 



w{zi) = 



WL{zi)r{zi) 



(8) 



The problem is that one needs to know the cross correlation 
parameter r to infer the true dark matter power spectrum. 
In the previous section we implicitly assumed r = 1, which 
is expected on linear scales. 

While we cannot simultaneously measure the stochas- 
ticity and bias, we can infer an integrated upper bound. 
Once we have solved for the weight function w, we define a 
residual SiS 



(9) 



In the absence of noise, this is the lensing power which was 
not predicted from the galaxies, and corresponds to the in- 
tegral of (1 — r{z)'^)w'j^{z) over the lensing weight. If this 
integral is small, it places an upper bound on the variation 
of r at any redshift. 

With this framework in mind, we can substitute the ac- 
tual parameters. We will use a linear power spectrum evolu- 
tion model, A'^{z) = A^(« = 0) D(zf /(1 + z)^. T he notation 
= k^P(k)/2n^ is taken from lPeacocS |l99i). 

The comoving angular diameter distance is 



dz 



where H(z) is the Hubble constant at redshift z: 



H{z) = Ho[{l + zf(Sl.mZ + 1) - Q.Az{z + 2)] 



1/2 



(10) 



(11) 



For the an gular diam eter distance x we use the fitting for- 
mula from lPenl Jl999ft . 

For dark matter the lensing weight is 



wl{z) = -Q.mHo g{z){l + z) 
where 

g{z) = x{z)J dzns(z) -^-^^ . 



(12) 



(13) 



ns{z) is the normalized distribution of source galaxies. The 
two dimensional dark matter map is then 

'ii0.,0y)= J S{dA{z)e^,dA{z)ey)wL(z)^dz. (14) 

Analogously, we define a weighted galaxy surface density 
with respect to a yet undetermined weight w{z) 



= / 5g{dA{z)6x,dAiz)6y, z)w(zi)dz. 



(15) 



In analogy to equation Q we solve for io(z) and infer the 
combination of bias and stochasticity given by equation Q. 



4 APPLICATIONS 

We consider parameters for the CFHT Legacy Survey 
JVan Waerbeke fc MeUie3l2003D . and its expected accuracy. 
We use a source density comparable to the VIRMOS- 
DESCARTES survey of 30 galaxies per square arcminutes 
with a median source redshift oi zq = 1. The mean ellipticity 
is taken to be e = 0.3. The survey area is 200 square degrees, 
which is half a percent of the sky. We neglect the Poisson 
noise from the cross correlation lens plane galaxies, which is 
much smaller than the lensing noise in each redshift slice. 

The LS lensing survey measures the projected two di- 
mensional dark matter distribution. In principle, the three 
dimensional power spectrum estimation requires inverting 
the Limber equation, which can be nu merically unstable . 
We will use the half weighted approach of iPen et alj i2003l) , 
which leads to mostly uncorrelated error bins in the three 
dimensional power spectrum. For simplicity, we use the 
power spectrum at the angular diameter distance of the 
typical source redshift zq — 0.4. The mapping becomes 
k = l/x{zo) ~ (i/lOOO)/!,"^ Mpc"\ We bin the power 
spectrum in factors of 2 in Z, and use Gaussian errors 
AC, = {Ci + C^)/^{2l + l)Mf,^y. 

The assumption of zero stochasticity is likely to break 
down when the galaxy clustering is non-linear. We model 
this conservatively as the stochasticity parameter being in- 
versely proportionate to the variance 



1 



1 + A2 



(16) 



Current data suggests that galax;ies and dark matter are 
probably better correlated than that jPen et alJ l2003l : 
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Figure 1. The expected lensing power spectrum for the CFH 
Legacy Survey. Error bars show the variance due to lensing. The 
upper and lower thin hnes show the variance if galaxy tomogra- 
phy is used. In our model, we assume stochasticity proportionate 
to non-linearity, which limits the tomography. Realistic stochas- 
ticities are smaller, so the procedure may work to smaller angular 
scales (larger 1). The dotted line is the lensing shot noise power 
spectrum. 

iHoekstra et alJl2002h . so tomography will likely work bet- 
ter than this estimate. The good news here is that one will 
be able to measure the stochasticity, and know the reliability 
of the inference. 

Figure Q shows the expected accuracy of the LS dark 
matter power spectrum measurement. To obtain the level 
of stochasticity in equation H16|l . we took the matter power 
spectrum ai z = 0.4. We assumed 1/4 radial modes along the 
line of sight. This is a rough estimate, based on a comov- 
ing radial size of 0.5h~^ Gpc for the reconstructed galaxy 
distribution. The two solid lines show the ac curacy with a 
matc hed redshift survey, for example CLAR iCarlson et alJ 
l20fl(t) . If only photometric redshifts are available, the im- 
provement can be limited by the square root of the number 
of radial redshift bins, if these are less than the number of 
radial modes. At / = 50, we need 12 radial bins to a redshift 
of 1, which may be achievable with photometric redshifts. 
With a lower fidelity, the errors will grow. 

There are other sources of errors which we have not 
included in the figure. There are inevitably errors on b{z, r). 
When stochasticity is small, one would in general expect 
the error on the bias to the smaller than sample variance. If 
the dark matter and galaxy maps differ by a multiplicative 
constant, one can just divide the maps and measure this 
constant perfectly, with no sampling variance. Measurement 
noise does limit at some level how well this can be done, 
which depends on both the lensing noise and the galaxy 
distribution noise. 

We see that we obtain significant improvement in the 
linear scale power spectrum estimation. The total statistical 
information is still dominated by the smaller non-linear scale 



power, for which many more modes can be measured. These 
non-Hnear scales in principle encode the power spectrum 
and its history, but require detailed N-body simulations to 
cafibrate. There are also questions of fundamental limits to 
the precision of non-linear scale power, since baryons back 
react on the dark matter. The distribution of baryons is 
intrinsically harder to predict since they can heat and cool 
in ways that is not predictable from first principles. 

Future lensing surveys in the next five years can cover 
the whole sky, and can map so urces at di stances beyond the 
epoch of reionization at z > 7 iPenll2003h . It is also possible 
to map the redshift distribution of galaxies to similar dis- 
tances on the whole sky with the Square Kilometer Array. 
This allows one to beat the normal sample variance limita- 
tion, also sometimes called 'cosmic variance'. In an interval 
of angular wavenumber / — 21 there are ~ 31^ two dimen- 
sional modes on the sky. This can in general be increased 
to ~ modes using tomography. But actual accuracy that 
can be gained will depend on the effects of non-linearity and 
stochasticity in the distribution of galaxies. 



5 CONCLUSIONS 

Gravitational lensing is a clean procedure to measure the 
power spectrum of dark matter at low redshifts. It is already 
providing competitive constraints on cosmological parame- 
ters, and does not make any assumptions about the nature 
of galaxy biasing nor any other gastrophysical process. 

Like the cosmic microwave background, gravitational 
lensing accuracy is limited by sample variance. It can only 
provide two dimensional power spectra. Some of the most 
important open cosmological questions involve precision 
measurements of the cosmic equation of state. For this pro- 
cedure to be successful, one needs accurate measurements of 
the dark matter power spectra, ideally in the linear regime 
and at different redshifts. This is a regime where gravita- 
tional lensing is significantly restricted in accuracy. 

We have shown in this paper how one can improve on 
the accuracy by using the cross correlation to galaxies with 
redshift information. One can reduce the errors without any 
assumptions on galaxy bias, even if it evolves and depends 
on scale. The only requirement is a low stochasticity. The 
total actual stochasticity of the data set is measurable, and 
can be folded into the error analysis. 

We estimated this improvement on two concrete sur- 
veys: the Canada-France-Hawaii- Telescope Legacy Survey 
Weak Lensing Survey, and the Canadian Large Adaptive 
Reflector deep survey. With the added distance information, 
one can reduce the error bars on linear scales Z ~ 50 by up 
to a factor of 3. 

I would like to thank Uros Seljak for helpful discussions. 
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